    % y = x .* log(sqrt((x.^2 - 1)) + x) ...
    %     - sqrt((x.^2 - 1)) ...
    %     - 0.5.*x;


% function [y] = fcode(x)
%     y = x.^2 .*sin(x.^2 - x - 2);
% end
% 
% 
% function [y] = ff(x)
%     y = -fcode(x);
% end
function [y] = fcode(x)
    y = x.^2 .* sin(x.^2 - x - 2);
end

% function y = fabcde(x)
%     y = 0;
% end